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In this article, we give a characterization of all the unitary transformations that can be synthesized 
in a given time for a two-qubit system in the presence of general time varying coupling tensor. This 
characterization helps to compute the minimum time and the shortest pulse sequence for generating 
a general two-qubit transformation under non-stationary interactions. The methods presented here 
can be applied in design of time optimal pulse sequences for transferring coherence and polarization 
between coupled spins with time varying couplings as in solid state NMR under magic angle spinning. 
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I. INTRODUCTION 

An important question in quantum information science 
is to determine the minimum time required to perform a 
quantum computation using a set of physical resources. 
Since two-qubit gates are the building blocks of quantum 
information processing, it is of fundamental interest to 
find the minimal time required to implement a unitary 
operation on a two-qubit system using the interaction 
Hamiltonian Hd, and the local unitary operations on the 
two qubits. This problem was studied in j3[, where it was 
shown that any two qubit unitary propogator Up can be 
expressed as 



Uf = U 2 (JJ V k cxp(-tH d t k ) Vl) U u 



(1) 



k=l 



where U\,U2,Vk are local unitary transformations and 
the effective Hamiltonians VkHdV^ all mutually com- 
mute. Under the assumption that the synthesis of local 
unitaries takes arbitrarily small time, the minimum time 
to produce a desired Uf is the smallest value of tj- 
in equation (JTJ) 3]. This characterization of time opti- 
mal trajectories is used in [6j to explicitly compute an 
elegant expression for the minimal time for synthesis of 
arbitrary unitary transformation of two qubits. Alter- 
nate proofs for time optimality have been presented in 
0, EH- There is now a considerable literature on the 
subject; see for example, ((1-IT4| and references therein. 

All these investigations assume that the interaction 
Hamiltonian Hd is fixed. In this paper, we consider 
the general problem when Hd varies with time. For ex- 
ample, in solid state NMR [l7|, the interaction between 
the spins are varying with time during magic angle spin- 
ning when the sample is rotated around an axis making 
an angle of 9 m = tan -1 (-\/2) with the static magnetic 
field B . As a result the dipolar couplings between nu- 
clear spins that have an orientational dependence of the 



form 3cos 2 (#) — 1 averages out( 9 is the angle of inter- 
nuclear axis with the static magnetic field), leading to 
better resolved NMR spectrum [17| ■ An important prob- 
lem in multi-dimensional solid state NMR experiments 
is to find radio- frequency pulse sequence that re-couple 
desired spins whose interactions are being modulated in 
time by magic angle spinning. Finding short pulse se- 
quences that transfer polarization or coherence between 
coupled nuclear spins under time varying interactions is 
of interest in solid state NMR. In this paper, we give 
a complete characterization of all the unitary transfor- 
mations that can be synthesized in a given time for a 
two-qubit system in presence of general time varying cou- 
pling tensor, assuming that the local unitary transforma- 
tion on two qubits can be performed arbitrarily fast (on a 
time scale governed by the strength of couplings). From 
the perspective of quantum control theory, this problem 
is equivalent to characterizing the reachable set of the 
Schrodinger equation 



U{t) = -i[H d {t)+Y j v j {t)H j ]U{t), 



(2) 



where U G SU(4) and Hd{t) is the interaction Hamil- 
tonian that is internal to the system and ^2j = xVj(t)Hj 
is the part of the Hamiltonian that can be externally 
changed, and generates the local unitary operations. 
We assume the control parameters Vj are a priori not 
bounded. 

Before stating the main result, we review some back- 
ground material. 



II. BACKGROUND 



A. Majorization 



For an element x = (x\, ...,Xk) T of R fc we denote by 



{x\, ...,x\.) T a permutation of x so that xj > if 



i < j, where 1 < i,j < k. 
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Definition 1 (majorization) A vector x £ R fe is ma- 



2 



jorized by a vector y G M. K (denoted x -< y), if 

E4<X>i w 



for fc = 1, . . . , D — 1, and the inequality holds with equal- 
ity when k = D. 

Proposition 1 x < y iff x lies in the convex hull of y 
and all its permutations Piy, where Pi are permutation 
matrices. 

Proposition 2 (Schur, Horn Q, [1]) For an element 
A = (Ai, A„) T , let D\ be a diagonal matrix with 
(Ai,...,A„) as its diagonal entries, let a — (a\, a n ) T 
be the diagonal entries of matrix A = K T D\K, where 
K G SO(n). Then a -< A. Conversely for any vector 
a -< A, there exists & K £ SO(n), such that (oi, a„) T 
are the diagonal entries of A — K T D\K 

Following 0-0], for an element x — (xi,x 2 , xs) T of R 3 , 
we introduce the vector x = (\x\\, \x 2 \, \x3\) T , and define 
the s-order version x s of x by setting x\ = x\, x 2 = x 2 , 
x 3 = sgn(x\x 2 X2,)x 3 . 

Definition 2 (0-@|) The vector x G R 3 is s-majorized 
by y G R 3 (denoted x ~< s y) if 

x\ < yt 

x s 1 +x s 2 +x s 3 <yi + y s 2 + yl (4) 

S t S S ^ S 1 s s 

x 1 +x 2 -x 3 <y 1 +y 2 - y 3 

B. Canonical Decomposition 

An arbitrary two-qubit Hamiltonian can be parame- 
terized 

H d (t) = Kg>(a(t)-a) + (b(t)-a)(g>I+^2M ij {t)a l (8>a 3 (5) 



where i,j G {x, y, z} and a = (a x ,a y ,a z ), b = (b x ,b y ,b z ) 
are real 3-vectors, M is a 3 by 3 real matrix, and 3 = 
{a x , a y , a z ) is the vector of Pauli operators. 
Let H' d (t) be the non-local part of Hd(t), i.e., 



H' d (t) 



i,3 



{t)a t <g> , 



Since we assume that the local unitaries can be generated 
in arbitrarily small time, all the unitaries transformations 
that can be synthesized in a given time under Hd(t) can 
also be synthesized under H' d (t) and vice versa Q. We 
therefore consider Hd(t) and H'Jt) are interchangeable 
resources under fast local unitaries. From now on we 
assume Hd(t) has only non-local terms. 



Proposition 3 (Canonical Decomposition [3;|, 0| ) 
Any two-qubit non-local Hamiltonian H can be written 
in the form 

H = {A®B)\e^ a x ®a x + 2 <j y ®<j y + 6 3 (j z ®a z ){A®B) 

(6) 

and any two-qubit unitary U G SU{A) may be written in 
the form 

U = {A 1 ® B 1 )e-^ e "' y ^^ +e "' y y®' J y +e " mz \A 2 ® B 2 ) 

(7) 

here A, A\ 1 A 2 , B, B\ 1 B 2 are singlc-qubit unitaries, and 



)?>e*>\ef\ 



(8) 



We call 9itj x <g> a x + 2 a y <g> a y + 6 3 cr z ® a z and 
e -i(eY**®*x+e%<Ty®<T v +ev*z®a z ) the canonica i f orm f H 

and U respectively, and 9 H and 6 17 the canonical pa- 
rameters of H and U respectively. For a 3-vector /?, we 
denote 

H0 = f3ia x ®a x + f3 2 <j y ®<r y + /3 3 a z ® a z 



u fj — a 



C. Magic Basis 

The magic basis is a vector space basis for two-qubit 
pure states: 

|00) + |11). .|00)-|11) 



V2 



V2 



J01) + |10) |01)-|10) 



V2 



V2 



(9) 



The basis change from the standard basis 
{|00), 1 01} , 1 10} , |11)} to the magic basis is given by 
<5 _1 , where 



Q =T2 



( 1 i 

i 1 

i -1 

V 1 -i 



For elements U G SU(4) the map U -> Q^UQ 
reflects the isomorphism between SU(2) (8 SU(2) and 
50(4) 16 ! j. When expressed in the magic basis, canon- 
ical form Hamiltonian and unitaries are diagonal. In 
magic basis, the canonical decomposition takes the form 
H d = K T D H K, U = RDuS, where K : R and S are 
real orthogonal matrices, and Dh, Djj are diagonal ma- 
trices. The diagonal elements of Dh and Djj are easily 
written in the terms of the canonical form parameters 
9\ie {H,U}. Define 



= V l +V 2 -V 3 ,(f 2 = »i 

4 = -oi+8i+oivi = - 



'3: 



(10) 
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The diagonal elements of Dh are {fi , tp%, (f^, <p±) 
and the diagonal elements of Djj are 
(e~ lVl , e~ lKP2 , e~ llf3 , e~ IV4 ). Under choice of magic 
basis, H d is real symmetric and {*-pf, (fi 2 , V>§ > <P±) are its 
eigenvalues. Eq.© together with Eq. pD]) implies that 

Proposition 4 ([7]) Let a and (3 be two real s-ordered 
3- vectors, let A and p be the 4- vectors related to a and 
j3 respectively via (|10|). then X ^. p iff a < s (3. 

The proof follows from the definitions. 

III. RESULT 

The main result of this paper is as follows: 

Theorem 1 Let H (t) be the canonical parameters of 
H d (t) in © and 9(T) = 6 H (t)dt, where the integra- 
tion is performed for each entry of the vector. All the 
unitary operators that can be generated within time T 
with Hd(t) and fast local unitaries are given by the set 

K(T) = {KiU$K 2 \Ki,K 2 e SU(2)®SU(2),0~< a 6{T)} 

Remark 1 We prove this theorem by using the choice of 
magic basis. In this basis, {—iHj} are skew-symmetric 
matrices and generate the group 5*0(4). The interac- 
tion Hamiltonian H d can be expressed as H d = K T D X K, 
where K € 50(4) and D x a diagonal matrix with diago- 
nal entry A related to 6 H via HTDJ}. Let <f(T) = J Q T X(t)dt, 
then in the magic basis 

K(T) = {Re- iD ?S\R, S £ 50(4), -< f (T)}. 

Proof: Under the choice of magic basis, we can write 
U(t) = R(t)A(t)S(t), where R(t),S(t) € 50(4) and A(t) 
be diagonal matrix. Assumption of fast local unitaries 
implies we can generate 50(4) instantly , so it suffices 
to prove all we can generate for the A part is e~ tD p , 
< l(T). 

Assume U(t) = R(t)A(t)S(t) is a trajectory of Eq.©, 
then A(t) = R T \t)U(t)S T (t) and A(t) = 

R T {t)U{t)S T {t) + R T (t)U(t)S T (t) + R T {t)U{t)S T {t). 

Let R T (t) = r(t)R T (t) and S T (t) = S T (t)s(t), substitut- 
ing for U(t), we get R T (t)U (t) S T {t) = 

m 

R T (t) [-iH d (t) - i^2v ] (t)H :j }R(t)(R T \t)U(t)S T \t)). 
i=i 

Using H d {t) = K T (t)D x K(t) we get R T (t)U (t) S T (t) = 

m 

R T (t) {-iK T (t)D x K(t) -iJ2^(t)Hj]R(t)(A(t)). 



Let P(t) = K(t)R(t) and denote h(t) = 
R T (t) [- i J2 7 ILi v j(t)H J }R(t). Equation for evolu- 
tion of A{t) then takes the form A(t) = 

r(t)A(i) + [P T (t)D x it)Pit)]A(t) + h(t)A(t) + A(t)s(t). 

(11) 

Notice that r(t),s(t), and h(t) are in so(4) (skew sym- 
metric matrices of dimension 4) and hence their diagonal 
entries are all zero. When multiplied by a diagonal ma- 
trix A(t), the diagonal entries remain zero. Therefore in 
the evolution equation of A(t), these terms must sum to 
zero and we can discard these terms. We get 

Ait) = D_ im Ait) (12) 

where pit) is the diagonal entries of P T it)D x ^P(t). 
Since we can generate elements of 50(4) in arbitrarily 
small time, R(t) and hence P(t) can take value of any el- 
ement in 50(4) and from Proposition pit) can take 
any element of the set {/T(t)|/J(t) -< A(t)}. 

From Eq.([I2]). we get A(T) = e~ iD $, where = 
J T p(t)dt, p{t) -< A(i). We first prove -< f(T), and 
then show that j3 can take on the values of any vector 
majorized by 7(T). 

E%-( T )=E [ T ^m (13) 



E^=E r^Wdi (14) 
3=1 3=1 Jo 

where tr is some permutations and k — 1,2,3,4. On 
subtracting Ea.(|I"4"]l from Ea. (|T3|) . we get 

E - E = f E 4 w - E ^o) (wis) 

3=1 3=1 - 70 3=1 3 = 1 

Since /x(t) x A(i), £* =1 A )(*) - E*=i /V(3) W > 0, and 
from Eq.®, £? =1 t£(T) - £*=i $ > 0. Obviously 
when k = 4, both terms equal 0, the equality holds, so 
/3 -< f(T). 

We now prove can take on the values of all the vectors 
majorized by 7(T) , which is the convex hull of 7(T) and 
all its permutations. If we take i2(t) = K T it), then f3 — 
7(T). It is also easy to see j3 can take all the permutations 
of 7(T), so we just need to prove that the vectors /3 can 
reach is a convex set. Let a E [0, 1], 

01= f piit)dt (16) 
Jo 

02= f ik(t)dt (17) 
Jo 
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then 

a0x + (1 - a)0i = [ afait) + {1 - a)p 2 (t)dt (18) 
Jo 

but api\{t) + (1 — a)p2(t) ~< A(i), so a/3i + (1 — a)/3i can 
also be achieved. Q.E.D 

Given these theorems, we can compute the minimum 
time needed to generate any unitary operator U in SU (4) 
with H<i(t) and fast local unitaries. 

Theorem 2 Using the Hamiltonian Hd{t) and fast local 
unitaries, a two-qubit gate U can be generated within 
time T iff there exists a vector n = (nx, 712, TI3) of integers, 
such that j3fi — 9 u + ?n satisfies 

/3ft / 9 H ^dt 
Jo 

where 9 u and 9 Hd ^ are the canonical parameters of U 
and Hd(t) respectively. The minimum time required to 
simulate U is given by the minimum value of T > such 
that either 

/3(o,o,o) -<* [ T 9 H ^dt (19) 
Jo 

or 

■<* f 9 H ^dt (20) 
Jo 

holds 

The proof follows the treatment in [(|. 

Proof: Recall that all commutators [aj ® <Jj,ak <8> 07-] 

vanish, and that exp(— if <Xj £§) <7j) = — i<Xj (S> Oj is a local 

gate. This implies that 9 u + %n represents all vectors 
compatible with the gate U. It is straightforward to check 
from Eq.fTjTJ) that for any two vectors x and y, with com- 
ponents xi > x 2 > \x 3 \, Vi>V2> \ys\, if Vi > 3xi, then 
x -< s y. By definition ^ > 9^ > 0, if some component rij 
of n fulfills \rij \ > 1, then the maximal component of the 
reordered version of 9 u + %n is at least which implies 
+ \n. Therefore we can restrict our attention 
to vectors n with \rij\ < 1. A case by case check shows 
that for n £ {(-1, -1, -1), (0,-1,0), (0,0, -1), (0,0, 1)}, 
9 U + f (-1,0,0) ^ s 9 u + |n, and for the remaining 

vectors n, 6 u + § (0, 0, 0) ~< s 6 u + ffi. Thus the result 
follows. 



IV. EXAMPLE 

We now work an explicit example on finding the 
minimum time to synthesize a desired unitary under 
time varying couplings. Assume the interaction Hd(t) 



takes the form D{t){a x ® a x + a y ® a y — 2a z ® a z ). 
We compute the minimum time to generate a swap 
gate corresponding to the unitary transformation U — 
exp—ij(a x ® a x + a y <8> a y + a z ® er z ). 

To fix ideas, consider the case when D(t) is constant, 
say D > 0. The canonical parameters of Hd(t) and !7 
are D(2, 1,— 1) and f(l, 1, 1) respectively. The mini- 
mum time to generate U is the minimum T that satisfies 
f(l,l,l) ^ s #7X2,1,-1) or f (1,1,1) + f (-1,0,0) ^ s 
DT(2, 1,-1), which is j^p- The strategy to generate U is 
to use selective excitation on first spin preparing an effec- 
tive Hamiltonian D(— a x ®a x + a y ®<J y J r2<j z ®<j z ), which 
evolves units of time. This is followed by evolution of 
effective Hamiltonians D(—2a x <E) cr x + <J y <E) <J y + a z ® a z ) 
and D{—a x ® a x + 2a y ® a y + <r z ® a z ) for units 

of time each. In the end we apply a local unitaries 

e -if<r x ®<Tx = _j CTx ^ Gx 

Now consider the time-dependent case, which models 
the variation of coupling strength between homo-nuclear 
spins under magic angle spinning (l7| . The dipolar inter- 
action strength D(t) during magic angle spinning varies 

in time as D{t) — J) 3cos ^f^ -1 , where 9{t) is the angle 
internuclear vector makes with the Bq field. The angle 
9(t) changes as the sample is being rotated around an 
axis making an angle 9m = tan -1 (v2) with the B$ field. 
Let (3 denote the angle internuclear axis makes with the 
magic angle axis. Then we can express 9(t) as 

cos(0(t)) = cos(/3) cos(#m) + sin(/3) cos(wi) sm(0M), 

where ui is the spinning frequency. D{t) is then a peri- 
odic function. We choose /3 = ? and plot modulation 
of D(t) in figure [TJ Each period of D(t) can be di- 
vided into two parts, {D(t) < 0}\J{D(t) > 0}. Let Si, 
S2 denote the area of these two parts respectively, i.e., 

Si = - I{D(t)< 0} D ^ Si = I{D(t)>0} D ^) AL We find 

that S 1 = S 2 = l^lD. 

The canonical parameters for Hd(t) are 

J D(i)(2,l,-1) for D(t) > 
\ -D(t)(2,l,l) for D(t) < 

i.e., (2|D(i)|, \D(t)\,-D(t)). Using theorem ©, we get 
the minimum time to generate U is the smallest T that 
satisfies 

^<j\\D(t)\-D(t)dt 

when w >> D, 3\D(t)\ — D(t)dt is approximately 
n(25i + 4S2), where n is the number of periods of D(t) 
within time T, so the minimum n = r 4 ( 2 gf+4s 2 ) 1 = 
|- 0.263 o) -| an( j ^.jjg mm i mum time T is approximately 
27r^. The pulse sequence prepares effective Hamiltoni- 
ans (— <t x ® CTj. + cT a ® aj, + 2a z ® cr z ), (— 2cr ;c ® o x + a y ® 
o-y + er 2 ® c z ) and (— ct^ ® cr^ + 2cr y ® a y + a z ® cr z ) 
for n/3 periods each, in the part of the period when 
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0.5 1 1.5 2 2.5 3 



Time (units of 1/ro) 

FIG. 1: The figure shows the modulation of the coupling 
strength D(t) as function of time for ft = ^ 

D{t) > 0. Similarly, we prepare effective Hamiltonians 

(— a x ®a x — a v ®a y +2a z ®(j z ), (a x ®(j x +2a y ®(jy—a z ®a z ) 



and (— 2a x ®a x + (Jy®(Ty + <J z ®<7 z ) for n/3 periods each, 
in the part of the period when D(t) < 0. As before, we 
apply a local rotation e~ l t^®^ in the end. 



V. CONCLUSION 

In this paper, we studied the problem of time-optimal 
synthesis of a unitary transformation for coupled qubits 
under non-stationary interactions. Under the assump- 
tion that local unitary transformations can be synthe- 
sized arbitrarily fast, we characterized the time optimal 
trajectories and the minimal time to prepare a general 
two qubit rotation under general time varying coupling 
tensor. These results generalize the results presented in 
[!) Hi Q for stationary coupling Hamiltonians to the non- 
stationary case. The problem considered in this paper 
was motivated by design of time optimal pulse sequences 
for controlling coupled spin dynamics in solid state NMR 
spectroscopy, where couplings between spins are modu- 
lated in time due to magic angle spinning. The results 
presented here are of fundamental interest and may find 
applications in some implementations of quantum infor- 
mation processing. 
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